Level density for deformations of the Gaussian orthogonal ensemble 
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Formulas are derived for the average level density of deformed, or transition, Gaussian orthogonal 
random matrix ensembles. After some general considerations about Gaussian ensembles we derive 
formulas for the average level density for (i) the transition from the Gaussian orthogonal ensemble 
(GOE) to the Poisson ensemble and (ii) the transition from the GOE to m GOEs. 
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I. INTRODUCTION 

Deformed random matrix ensembles were introduced 
by Rosenzweig and Porter 0, Q to classify the conser- 
vation of electronic spin and orbital angular momentum 
in the spectra of complex atoms. Many varieties of de- 
formed ensembles have since been constructed. Recent 
reviews of deformed ensembles (they are also called tran- 
sition ensembles) can be found in Refs. 0, E|- Earlier 
results on deformed ensembles were reviewed in Ref. [fj. 
Amongst the applications they have found we mention 
the breaking of time reversal invariance and the break- 
ing of symmetries . 

In this paper we are concerned with the level density 
for deformed Gaussian orthogonal random matrix ensem- 
bles. Although the level density of the standard random 
matrix ensembles is not directly related to the physical 
many particle level density it is essential to the proper un- 
folding of fluctuation measures . Unfolding is a trans- 
formation which leaves sequences of energy levels with a 
constant average density. Normally, experimental energy 
levels and the energy levels obtained from random ma- 
trix theory are unfolded (which removes the secular en- 
ergy dependences) before fluctuation measures for both 
are calculated and compared. While the level densities 
of most random matrix ensembles are different, the be- 
haviour of the fluctuation measures of a wide class of 
ensembles is the same and hence denominated universal. 
For instance, each of the non-Gaussian ensembles stud- 
ied in Ref. [!j is shown to have a characteristic average 
level density, however, after unfolding, the nearest neigh- 
bour spacing distributions and number variances of all 
are found to be identical to those of the standard Gaus- 
sian random matrix ensembles fi"o| . 

In the following we derive formulas for the level density 
of two deformations of the Gaussian orthogonal ensem- 
ble (GOE). The first describes a transition from the GOE 
to an ensemble with Poisson fluctuation statistics, that 
is, from the GOE to an ensemble of diagonal matrices 
whose elements are independently Gaussian distributed 
[ill Il2| . Such random matrix models are of interest to 
the analysis of experimental data |13( and of dynamical 
models [l4| whose fluctuation properties are intermedi- 
ate between Poisson and GOE. The second deformation 
describes the transition from the GOE to m GOEs, that 
is, from the GOE to a block dia gon al matrix with m 
blocks each of which is a GOE [Hlll[l3- If the blocks 



are labelled by quantum numbers such as angular mo- 
mentum pj or isospin |18| then the transition ensemble 
classifies their conservation or non-conservation. The m 
GOE to GOE transition is also relevant to the analysis 
of symmetry breaking in quartz blocks L9l l2Ct |21| . 

Our method extends that of Wigner |22} who showed, 
by assuming that terms containing patterns of unlinked 
binary associations dominate the averages of the traces 
of powers of matrices, that the level density for certain 
random matrix ensembles could be simply expressed as 
a Fourier transform (see also Sec. Ill D Ref. Previ- 
ous results for the level density of deformed GOEs have 
been derived using Stieltjes transform methods |2;i 1 12 11] . 
These methods are in fact general enough to treat de- 
formations of, and interpolations between, any class of 
matrix ensemble. However, the formulas obtained in the 
present paper have the advantage that that they are ex- 
plicit and simple to evaluate numerically. Other methods 
for calculating the level density (2^, |26| are restricted to 
deformations of the Gaussian unitary ensemble (GUE). 



II. INTERPOLATING GAUSSIAN ENSEMBLES 

The joint probability distribution of the matrix ele- 
ments of a matrix, H, for an interpolating Gaussian en- 
semble may be expressed as 

P{H,A,B) = Z- 1 (A,B)exp [- (AtrH 2 + Btr Hi)} , 

(1) 

where Z is a normalization factor and tr H denotes the 
trace of H . The structure of the matrix H\ is chosen 
in such a way that it defines a subspace of H . The pa- 
rameters A and B define the (energy) scale and degree 
of deformation. When B — > the joint distribution be- 
comes 



P(H,A,0) = Z- 1 (A,0)exp(-AtriJ 2 



(2) 



If we assume that H is real symmetric then the variances 
are given by 



« 4A 



(3) 



so that the limit B — > defines the GOE. When B — > oo, 
the elements of H\ vanish and H is projected onto a 
sparse matrix, Hq, whose elements are on the comple- 
mentary subspace of H±. Therefore, the random matrices 



2 



generated by Eq. JI} are the sum of two terms 

H = H + H 1 , (4) 

with the variances of the matrix elements of Hq given by 
the r.h.s. of Eq. © and those of Hi given by 



H 2 _ t 1 
1%3 4(A + B)' 



(5) 



This shows that when B goes from zero to infinity, the 
ensemble undergoes a transition from GOE to the Gaus- 
sian ensemble of sparse matrices defined by the choice 
of the structure of Hi or, equivalently, of H . It is also 
instructive to introduce the parameter 



a = (l + B/A) 



-1/2 



(6) 



which measures the relative strength of Hi and Ho- The 
transition from a—0 to a— I corresponds to the transition 
from B=oo to B=0. 



III. LEVEL DENSITY FOR THE GOE 

We now proceed to give a derivation of the semicircle 
law, valid for the GOE, which is equivalent to Wigner's 
for the ensemble of random sign symmetric matrices |22j . 
In general, the average level density (ALD) may be writ- 
ten as a Fourier transform, 



where 



/oo 
dEp(E) 
-co 

y> (ik) n 



ih-E 



dEp(E)E n . 



(7) 



(8) 



(9) 



Let Ek denote the eigenvalues of an N x N matrix H, 
which satisfies Eqs. J2J and ©■ From the exact expres- 
sion, 



i N 

(E) = ^6(E-E k ), 



(10) 



fc=i 



for the ALD, one obtains the following connection be- 
tween the moments of the eigenvalue distribution and 
the moments of the matrix elements: 



E h 



1 



dEp (E) E n = ^H r ' 



Substituting Eq. lfTT)l into Eq. JSJ we obtain 



N ^— ' n! 

n=0 



(11) 



(12) 



It is also useful to note that by differentiating Eq. (JSJ) the 
moments of the eigenvalue distribution can be expressed 
as derivatives of the Fourier transform of the ALD eval- 
uated at zero, that is, 



(13) 



where fW(0) = [d n F/dk n ] k=Q . 

For large matrices, N ^> 1, the average of the trace 
is dominated by the terms in which matrix indices may 
be contracted in such a way that we are left with s = ^ 
pairs of matrix elements. This means that we can write 



txH n ~ C s N s+l ct Zs , 
where we have introduced the notation 

^ 1 

= 4l 



(14) 



(15) 



for the variance of the off-diagonal matrix elements. The 
factor C s counts the number of contractions that leads to 
pairs. We observe that in a contraction s — 1 indices are 
eliminated and s + 1 remain which explains the power of 
N in the above expression and allows us to write 



(2s)! 



s(s-l)!(s + l)! : 



(16) 



where the binomial factor counts the number of ways s+1 
indices can be extracted out of 2s ones. The factor is then 
divided by s, that is, the number of ways s — 1 indices 
can be eliminated without changing the contraction. 

Substituting Eq. O m Eq- dPi we find that the 
Fourier transform of the average density is given by 



Fi (a, k) 



1 X 



(-Y(ka/2) 
s!(s + l)! 

Ji (ka) , 



2s+l 



ka/2 



(17) 



where Ji (x) is a first-order Bessel function and 

a = y/N/A. (18) 

Using the formula 



(l-x 2 ) 1/2 = 



y-'Ji^e-^dy, 



(19) 



the Fourier transform in Eq. Q may evaluated to obtain 
the semi-circle law, 



Pi (a,E) 



0, 



E > a 



(20) 



where the radius of the semi-circle is given by Eq. JTSJl. 
The cumulative level density (CLD), x(E), is defined by 



x(E) = / p(E')dE', 
Jo 



(21) 
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FIG. 1: Graph of p a (E) for several values of a with iV=1000. 
The calculations where performed with A=N/A for which the 
radius of the Wigner semi-circle is a=2. The solid lines were 
calculated using Eq. I|36|l and the dotted histograms by nu- 
merically diagonalizing an ensemble of 100 matrices. 



which for the semi-circle law is found to be 

r -1/2, E < a 

X 1 (d,E)=l EVa^-E^+a^Tcsin(E/a) _ |£| < fl ) (22) 

1 1/2, £ > a 

The average number of levels up to energy E is given in 
terms of the CLD by N[x(E) - sc(-oo)]. 



IV. THE TRANSITION GOE TO POISSON 



where E j and \ j) are defined by the eigenvalue equation 

— Eoj\j) and Eq denotes the ith moment of the 
eigenvalue distribution of Hq. Eqs. (j2"3*|) and l (2*5)l allow 
us to write [cf. Eq. ijTT])] 

2 /' OC /» CO 

-tr~ff^« / dEopo(Eo) dEipi(Ei) (E + E x ) n , 

•/ — OO </ — OO 



(26) 



where po and pi are the average level densities corre- 
sponding to i/o and ifi respectively. Eq. (|2l))l is only 
approximate because in resumming the binomial series 
we have kept linked as well as unlinked binary associa- 
tions of Hq and H\ . In the following subsection we show 
that for the case of the transition GOE-Poisson the re- 
sulting discrepancy in the ALD is small. We also show 
how to correct for this discrepancy. 
Substituting $Zj$ into jQ , wc find 

dEoPo(Eo) / dE lPl {E 1 )e^ Ea+E ^ 

-oo J —CO 

= FoWF^k), (27) 

where Fq and F\ are the Fourier transforms of pq and 
pi respectively. The ALD is obtained by substituting 
Eq. (|27|l into Eq. Q . Another representation is obtained 
by noting that since Eq. (|27|l is a product of Fourier trans- 
forms, the ALD of H is given by the convolution of the 
average level densities of Hq and Hi : 



P(E) 



dE'p (E')pi(E - E') 



(28) 



The only assumption required to derive Eqs. (|27|) and 
(|28|l is that the matrix elements of Hq and Hi be statis- 
tically independent. 



Before considering the transition from Poisson to GOE 
we derive a formula for the ALD for a more general case. 



A. Level density for Ho + Hi 

The trace of the nth power of Eq. (J1J can be written 



as 



trH n = tr(H +Hi) n 



E 



11 {n-l) 



tr H l H?~ l . 



(23) 



If Hq and Hi and are statistically independent we can 
write 



tvH l Q Hr l = J2 E Oj(i\ H i~ l \3) (24) 

(25) 



E l trHr l : 



B. Level Density for the transition GOE to Poisson 

To specialize the results of the last subsection to the 
transition GOE-Poisson, Hq is chosen to be the diagonal 
matrix Hq^ = EQiSij whose eigenvalues, Ecu, are inde- 
pendent random variables with Gaussian distribution 



Po(E) = 

The variance of Hq is thus 

5 



H 2 



2A 



2N 



(29) 



(30) 



We choose Hi to be a diagonal-less matrix whose matrix 
elements are independent Gaussian variables with zero 
mean and variances 



H 2 



,1 



AA 



i i 
a a 



1 



AN 



(31) 
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The ALD for H alone is given by Eq. (|2*9"|l and the 
CLD, Eq. PI)), bv 

x (E) = ~erf(Vl£). (32) 

The Fourier transform of Eq. I|29[l is 

F (k) = e- k2 /± A . (33) 
Note that for fixed a,^4-^ooasiV— >oo, so that 



Po(E) 



N- 



5(E). 



(34) 



The ALD for Hi alone is given by pi(aa, E), that is, by 
the semicircle law [see Eq. ()20|l] with the radius modified 
(a — * aa) in accordance with Eq. (|3 1 1> for the variance 
of Hi [cf. Eq. J3J]. It's Fourier transform is Fi(aa,k) 
[see Eq. (|f 71 ]. The ALD for Hi is given by the semicircle 
law in spite of the missing diagonal because off-diagonal 
matrix elements dominate (for N large enough). 

Using Eqs. and J7J) the ALD which interpolates 
between the Gaussian and the semicircle as a varies be- 
tween and f is then found to be 



2tt 



dkFo(k)Fi(aa, k)e 



-ikE 



7T A ./n X A 



(36) 



where we have changed the integration variable to x = 
aka and introduced the parameter 



A = a\/N. 



(37) 



In Fig. we compare Eq. (|31)|) for p a (E) with histograms 
constructed by numerically diagonalizing random matri- 
ces for several values of a (see figure caption for details). 
It is seen that good agreement is obtained, there be- 
ing however a small discrepancy in the transition region 
which is especially apparent in the graph for a=0.05. The 
CLD, Eq. J2U, may be expressed using Eq. lj36*|) as 



x a (E) = - /°°^ e - 2 /4A 2 Ji(x)sin 

7T Jo X 



AxE 



A 



(38) 



This equation provides a more accurate manner of un- 
folding fluctuation measures than the polynomial unfold- 
ing used in Ref . |27j • Eq. (|38|l is plotted for several values 
of a in Fig. [Z\ 

The ALD may be expressed alternatively using the 
convolution formula, Eq. (|28(l . In particular, when N — > 
oo for fixed a and a we find from Eq. and Eq. 
that 



Pa(E) 



N- 



pi(aa,E). 



(39) 




FIG. 2: Graph of x a (E) calculated using Eq. I|38|l for several 
values of a with iV=1000. The calculations where performed 
with A=N/A for which the radius of the Wigner semi-circle 
is a— 2. 



For finite N Eq. $Zj$ yields 



(35) Pa ( E ) = 



2VN 



aa+E 



aa+E 



I 3- ^a 2 a 2 -(E- E') 2 dE'. 

It it clear from these explicit formulae that the tran- 
sition parameter is a 2 N as was already argued in the 
original paper of Rosenzweig and Porter [l| . The model 
described in this section can be cast in the form 



H = Hn 



A 



Hi, 



(41) 



with c = 1/2, by modifying the definitions of the vari- 
ances, Eqs. I|31|) . Ref. [26j considered the statistics of 
ensembles of the form of Eq. (|41|l for arbitrary c. In par- 
ticular the asymptotic behavior of the ALD as N — > oo 
depends critically on c, Eq. (|39|l being valid only for the 
special case c = 1/2. 

Eq. (|36|l may be improved by comparing the resulting 
lowest moments of the eigenvalue distribution with the 
average trace of the corresponding powers of the Hamil- 
tonian. Considering the latter first we have 



tr H 2 



trH 2 



trHf 



NEq + CiN 



2«V 

AN 



(42) 



and 



tri? 4 = 



trff 4 



4tr H%Hl 
NEj + AEjCiN 2 



trHf 

2^2 



or a 
AN 



C 2 N 3 



16N 2 



N 
8A 2 



(A 4 + 4A 2 + 6) . 



(43) 
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FIG. 3: Graph of p a {E) showing the improvement obtained 
by demanding that the second and fourth moments of the 
eigenvalue distribution be exact for a=0.05 and iV=1000. The 
solid lines result from using Eq. 1361 and the dot-dashed lines 
from using the same equation modified in accordance with 
Eqs. (1481 and 1491 . The dotted histograms were obtained by 
numerically diagonalizing an ensemble of 100 matrices. 



FIG. 4: Graph of the peak values of the density p a (0), 
Eq. 13611 . as a function of the transition parameter A, 
Eq. $7}. For comparison we also show the limiting cases 
po(0), Eq. 1291 . and pi (era, 0), Eq. 1391 . as well as an interpo- 
lation formula given by Persson and Aberg |2sT] ol A . Ea lETty . 
Finally, we show p a (0) calculated using Eq. 1361 modified in 
accordance with Eqs. 1481 and 1491 . 



Inserting Eq. (|36|l into Eq. for the moments of the 
eigenvalue distribution we find that 



IF 1 



8A> 2 



2), 

(A' 4 + 6A' 2 + 6) 



(44) 
(45) 



Note that the coefficient of A' 2 in Eq. (|43|l is 6. This is a 
result of the use of Eq. I|26|) which included linked binary 
associations of Hq and H%. Eq. demands that 



£ 2 = _L trjff2; 



E 4 = —trH 4 . 

N 



(46) 
(47) 



Solving Eqs. igBJ) and (03) for A' and A' we see that 
Eq. will give the 2nd and 4th moments of the eigen- 
values distribution correctly if we make the substitution. 



1 
A 



1 _ 1 

A' ~ A 

A' 2 = A 2 



1 + A 2 



1 - (1 



4 
A 2 



A 2 



1 + A 2 



1 



(48) 



(49) 



In Fig. |3 we show the improvement to the ALD obtained 
by using Eqs. I|48|l and l|49|l in Eq. (|36|l for the worst case 
of Fig. n (a=0.05). We see that the modified formula 
gives the ALD essentially exactly. 



In Fig. 0] we show the peak value of the density of 
states as a function of A (see figure caption) . We see that 
Eq. (|36|l in both its modified and unmodified versions 
agree at the A=0 limit. We also see that both versions 
of Eq. (SHI) reach the N — > oo limit [Eq. (JSUJ] at roughly 
the same value of A (A«6). However the two versions de- 
viate from each other in the transition region, the largest 
difference occuring around A«1.5 (a«0.05 for iV=1000). 
Also shown for comparison is the interpolation formula 
for the ALD of Persson and Aberg : 



N2 



4aNi + 7N- 



-1.5a 



(50) 



THE TRANSITION GOE TO m GOES 



We now calculate the ALD for the transition from the 
GOE to a superposition of m GOEs. To proceed we 
again consider an ensemble of matrices of the form of 
Eq. |0J. Now H is a block diagonal matrix consisting 
of m blocks whose dimensions are Mi, i — l,2,...,m, 

m 

with Mi = N. The elements of Hq have zero mean 

i=l 

and variances given by Eq. ©. We define Hi to be zero 
where H is non-zero and elsewhere its elements have 
zero mean and variances given by Eq. JSJ. 
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Although the authors were unable to obtain an equation 
analogous to Eq. for higher powers of H , Eqs. 
and strongly suggest that 

»=i 

for large N. 

Substituting Eq l|57|l into Eq. (|12fl and again using 
Eqs. (jnjl, O and JTSJl we obtain 

m M- 

P^E^Px (58) 

i=l 

where pi is given by Eq. (|20() and 



FIG. 5: Graph of p a (E) illustrating the transition 3 GOEs 
to GOE. The three block have size M 1 =50, M 2 =400 and 
M 3 =550. The densities for a= 0, 0.7, 1.0 are shown. The 
calculations where performed with A=N/4 for which the ra- 
dius of the Wigner semi-circle is a=2. The solid lines were 
calculated using Eq. 158H and the dotted histograms by nu- 
merical diagonalizing an ensemble of 100 matrices. 



To obtain a formula for the ALD we note that 

JV 



trlP = ]T H% 

j,k=l 

N * 



with 



AA 



i=l 



— + a 2 (l -) 

N V N' 



(51) 
(52) 

(53) 



Considering a single line of H the of consist of a term 
which is the product of the variance of a single non-zero 
element of Hq with the probability for being in block i 
plus a term which is the product of the variance of a 
single non-zero element of Hi with the probability for 
being outside block i. In the sum in Eq. Il52[l the of are 
weighted with the fraction of lines which find themselves 
in block i. 

For the fourth power of H we find 



JV 



tvH 4 



2 E H % H l 

j,k,i=i 



N 

= 2 £ 

7=1 



fc=i 



2N 3 Y^at. 

i=l 



(54) 

(55) 
(56) 



1 

A 



M M 
N y N ' 

M l + X\l-f) 



(59) 
(60) 



From Eq. (JSBJ it follows that the CLD, Eq. (J2J, is given 
by 



N 



(61) 



with a:i given by Eq. (|22|) . 

In Fig. we compare Eq. I|58|l for the ALD with nu- 
merical simulations for several values of the a (see figure 
caption for details) and excellent agreement is obtained. 



VI. CONCLUSION 

In conclusion, by assuming that terms containing pat- 
terns of unlinked binary associations dominate the the 
averages of traces of powers of matrices, we have derived 
formulas for the average level density for two deforma- 
tions of the Gaussian orthogonal ensemble. The first de- 
scribes the transition from the Gaussian orthogonal en- 
semble to the Poisson ensemble and the second the tran- 
sition from the GOE to m GOEs. The formula obtained 
are in excellent agreement with numerical simulations. 
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